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We demonstrate that the ground state of a trapped spin-1 and spin-2 spinor ferromagnetic Bose- 
Einstein condensate (BEC) can be well approximated by a single decoupled Gross-Pitaevskii (GP) 
equation. Useful analytic models for the ground-state densities of ferromagnetic BECs are obtained 
from the Thomas-Fermi approximation (TFA) to this decoupled equation. Similarly, for the ground 
states of spin-1 anti-ferromagnetic and spin-2 anti-ferromagnetic and cyclic BECs, some of the spin 
component densities are zero which reduces the coupled GP equation to a simple reduced form. 
Analytic models for ground state densities are also obtained for anti-ferromagnetic and cyclic BECs 
from the TFA to the respective reduced GP equations. The analytic densities are illustrated and 
compared with the full numerical solution of the GP equation with realistic experimental parameters. 
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I. INTRODUCTION 

The advent of the optical traps paved the way for 
the first realization of a Bose-Einstein condensate (BEC) 
with internal spin degrees of freedom [T], also known as 
a spinor BEC. Since then, a lot of theoretical and ex¬ 
perimental studies have been performed on the spinor 
BECs [2Hi]. In contrast to a scalar BEC, which is char¬ 
acterized by a single interaction parameter, the spin- 
1 [5J and spin-2 [SJ '7\ BECs have, respectively, two 
and three interaction parameters. Depending on the 
relative strength of the interaction parameters, a spin- 
1 BEC in the absence of external magnetic field can 
be either in a ferromagnetic or an anti-ferromagnetic 
phase [5]. In the presence of magnetic field the ground 
state phase diagram of spin-1 condensate has been in¬ 
vestigated both for uniform [EHEI and trapped systems 
[IQ] HI. ■ Similarly, in the absence of external magnetic 
field, a spin-2 BECs can be in one of the possible three 
ground state phases: ferromagnetic, anti-ferromagnetic 
and cyclic 0 U\- The spin-1 and spin-2 BECs are de¬ 
scribed by three- and five-component complex order pa¬ 
rameters, respectively, thus leading to coupled mean- 
field Gross-Pitaevskii (GP) equations involving three- 
and five-component wave functions, which, unlike in a 
scalar BEC, could be complex in general. A numerical 
solution of these equations could be cumbersome for both 
spin- 1 HM3 and spin-2 mug BECs. In this paper, we 
propose simple and useful analytic models for the den¬ 
sities of ground states of quasi-one-dimensional (quasi- 
ID), circularly-symmetric quasi-two-dimensional (quasi- 
2D) [IB], and spherically-symmetric three-dimensional 
(3D) spin-1 and spin-2 spinor BECs. Here, we 
consider nearly-overlapping spatially-symmetric ground 
states only. Phase-separated spatially-asymmetric pro- 
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files [T]J] do not appear as ground states and will not be 
considered. 

The two interaction parameters for a spin-1 BEC are 
Co oc (ao + 2 d 2)/3 and ci oc (<22 — cio)/3 [5], whereas 
the three interaction parameters for a spin-2 BEC are 
Co oc ( 4 d 2 + 3 d 4 )/ 7 , Ci oc (04 — 0 - 2 ) 17 , and C 2 oc (7ao — 
10a 2 +3a 4 )/7 [H[7] , where ao, a 2 , and 04 are s-wave scat¬ 
tering lengths in total spin f tot = 0, 2, and 4 channels, 
respectively. For a ferromagnetic BEC, e.g. for Ci < 0 
for a spin-1 BEC, and for Ci < 0 and c 2 > 20ci for 
a spin-2 BEC, we find that to a very good approxima¬ 
tion the densities for different spin components irif of 
the ground-state wave function with magnetization m are 
multiples of each other. This allows one to replace the 
coupled GP equation for the ground-state wave function 
by a single partial differential equation, which we call the 
decoupled-mode (DM) equation. On the other hand, for 
an anti-ferromagnetic BEC, e.g. for Ci > 0 for a spin-1 
BEC, and for c 2 < 0 and c 2 < 20ci for a spin-2 BEC, we 
find that the densities, for some of the spin components, 
of the ground-state wave function with magnetization m 
are identically zero, thus reducing the original coupled 
GP equation to a system of two coupled equations for 
any non-zero magnetization. Similarly, for a cyclic BEC, 
e.g. for ci > 0 and c 2 > 0 for a spin-2 BEC, the five- 
component GP equation reduces to a system of two or 
three coupled equations. These reduced GP equations 
and the DM equation for the ground state of a spinor 
BEC in different parameter domains, valid in all spa¬ 
tial dimensions, are solved in the Thomas-Fermi approx¬ 
imation (TFA) (or local-density approximation) to yield 
simple analytic models for the ground-state densities of 
spin -1 and spin-2 spinor BECs. 

The TFA is applicable when the interaction energy in 
the GP equation is much larger than the kinetic energy 
term, so that the latter could be neglected, thus leading 
to simple analytic formulae for the condensate densities 
m■ In a repulsive scalar BEC, applicability of TFA 
requires that the size of the condensate R is much larger 
than the oscillator length Zq, he. R/Iq >> 1 [TS][TS]- The 
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spatial extent of the BEC in units of Iq is 


II. REDUCED MEAN-FIELD EQUATIONS 


_ R / Na\ 1/(:D+2) 

' h ~ \~v) 


(1) 


where V = 1,2,3 is the dimensionality of the space 
ng HU- The criterion is satisfied if the dimensionless 
parameter Na/lo >> 1. The ratio Na/lo is a measure of 
the strength of repulsive interaction. For a spinor BEC, 
the applicability of TFA for m / component requires that 
its spatial extent is much larger than Iq. 

There have been few studies to include the neglected 
kinetic energy contribution in the TFA CHI EH. Earlier, 
the TFA was used to study the ground state properties of 
binary condensates [2T1 and spin dynamics in quasi-ID 
spin-1 condensate [22J. Spin-orbit-coupled pseudospin- 
1/2 BECs under rotation have also been theoretically 
investigated using the TFA j23j . 

We use the experimentally realizable trapping poten¬ 
tial and interaction parameters to illustrate the present 
analytic models for ground-state densities in different pa¬ 
rameter domains. In the case of a spin-1 BEC, the back¬ 
ground scattering lengths of 8 7 Rb and 23 Na fall in the fer¬ 
romagnetic [23] [25] and anti-ferromagnetic [26] domains, 
respectively, and we use these to study the ground state 
properties. In the case of a spin-2 BEC, we employ 23 Na 
and 83 Rb BECs for the illustration. The background 
scattering lengths of spin-2 23 Na and 83 Rb correspond 
to the anti-ferromagnetic and ferromagnetic phases, re¬ 
spectively [6j. By tuning one of the scattering lengths 
of 23 Na, one can move from anti-ferromagnetic to either 
ferromagnetic or cyclic phase. Experimentally, such a 
change of single scattering length can be achieved by 
exploiting magnetic m and optical [28 Feshbach res¬ 
onance techniques. The results of the analytic models 
are also validated by a numerical solution of the origi¬ 
nal mean-field GP equations for quasi-ID and quasi-2D 
traps. 

In Sec. [IT] we present the full mean-field GP equations 
for spin-1 and spin-2 BECs and derive the reduced mean- 
field GP equations for the ground-state wave function in 
the different parameter domains. By assuming that the 
component wavefunctions of a ferromagnetic BEC are 
proportional to each other, which is indeed the case as 
suggested by numerical simulations, we derive the DM 
equation. By minimizing the ci- and C 2 -dependent en¬ 
ergy terms for the ground-state wave function, we obtain 
the reduced GP equations in all parameter domains. In 
Secs. |III| and [TV| we obtain the analytic models for spin-1 
and spin-2 ground-state BECs, respectively, by employ¬ 
ing the TFA to the reduced GP and the DM equations. 
A comparison of the analytic densities with the numer¬ 
ical densities obtained from the full GP equations leads 
to a very satisfactory agreement. In Sec. [V] we present a 
summary and concluding remarks. Some of the technical 
details about the derivation of the DM equation and the 
reduced GP equations in different parameter domains are 
presented in Appendix A and B. 


A. Spin-1 BEC 


The coupled GP equations for different spin compo¬ 
nents rrif = ±1,0, for a spin-1 BEC of N atoms of mass 
M each can be written in dimensionless form as [2] 


M±i0±i(x) = Hcp± i(x) + c 0/ 9</±i(x) ± c 1 F z <j)±i(x) 

± ( Cl /V2)F T 0 o (x), (2) 

Mo</>o(x) = H0o(x) ±c 0 p</o(x) + (ci/\/2)[i r _0_ 1 (x) 

+ F 1 + </)-|_i(x)], (3) 


where 


F ± =F X ±F y = V / 2[(/>± 1 (x)^ 0 (x) + c/>o(x)</> T i(x)], 

1. 


f z =p+i(x) — p_i(x), n = 


- -V 2 
2 


-U(x) 


(4) 

(5) 


and the component density pj = |</>y (x) | 2 with j = ±1,0, 
the total density p = JT p.j. and p±i, Po are the respec¬ 
tive chemical potentials and * denotes complex conju¬ 
gate. In 3D, the interaction parameters, Laplacian, and 
trapping potential are defined as 


co =- 


47r./V(ao ± 2d2) 
^0 ’ 


Cl = 


47rA^(a 2 - a 0 ) 

3Zn 


( 6 ) 


72 = &_ + + v( , = s 2 + ffV± 7 2 * 2 (7) 

dx 2 dy 2 dz 2 ’ [ ’ 2 ’ 1 J 


with x = {x,y, zj. Here l 0 = y/h/(Mw x ), = u y /ui x , 
7 = oj z /u x , where Lu x ,u] y ,uj z are the confining trap fre¬ 
quencies in x,y, z directions, respectively. When the 
trapping frequency along one axis, say co z , is much larger 
than the geometric mean of other two, i.e., u) z yjw x w y , 
then one can approximate the Eqs. © -(© by quasi 
two-dimensional (2D) equations which can be obtained 
by substituting m 


2 N^2Tr(a 0 ± 2 a 2 ) 2N\Z27r(a 2 - a 0 ) 

Co= - W, - >Cl = -’ (8) 

2 d 2 d 2 x 2 + P 2 y 2 

v + y(x) = —2-(9) 

in Eqs. © -©, here l z = \Jfif (Mui z ). Similarly, if the 
trapping frequencies along two axes, say y and z, are 
much larger than the third frequency oj x , Eqs. © -© 
can be approximated by quasi-ID equations which can 
be obtained by substituting 


co = 


2 N(ao ± 2a 2 )lo 


Wy Z 


Cl = 


2N(a 2 — ao)lo 


3 P yz 


f)2 2 

V2 = 5^’ n*) = y.x S *, 


( 10 ) 

( 11 ) 


where l yz = y/h/(Mui yz ) and ui yz = y/uj y oj z . Here length 
is measured in units of Iq, density in units of Iq V and 
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chemical potential in units of Hcu x , where V = 1, 2, 3 is the 
dimensionality of space. The total density is normalized 
to unity f dxp(x) = 1. The volume element dx = 2dX in 
ID, 2irXdX in 2D with circular symmetry, and 4TvX 2 dX 
in 3D with spherical symmetry, where X = |x| is the 
length of the vector x. In this paper, we consider isotropic 
3D and isotropic quasi-2D traps, i.e., /3 = 7 = 1 for 3D 
traps and /3 = 1 <C 7 for quasi-2D traps. 

Numerical calculation for the ground-state densities of 
a of a ferromagnetic BEC (ci < 0) has revealed that 
the component densities are essentially multiples of each 
other. This opens the possibility of writing a single de¬ 
coupled mode (DM) equation for the wave-function cj )dm 
for the ferromagnetic BEC and obtain the component 
wave functions as multiples of this wave function accord¬ 
ing to 


<M X ) = ay 0 dm (x), j = ± 1 , 0 , 


( 12 ) 


where ads, in general, are complex numbers. The con¬ 
ditions (12 1 when substituted in Eqs. Q and ^ lead 
to three different equations for the same wave function 
cf> dm - A consistency requirement on these three equa¬ 
tions leads to the single decoupled-mode (DM) equation 
for the wave function </>dm: 


^dm(x) = 


- -V 2 + E(x) +C^m( x ) 


^dmW, (13) 


with C = Ci = Cq + Ci and normalization f = 

1 , provided that 


1 ± m 


l“±i| = 


l“ol = 


\/l — n 


= j dx[p +1 - p_ 1 ] = |a+i | 2 - |a_i| 


(14) 

(15) 


where to is the magnetization. An equation similar to Eq. 
(13 1 with C = Co ~ (ao + 202 )/3, known as the single¬ 
mode approximation (SMA) [2H], was obtained before as 
an approximation to Eqs. ( 0 i- 0 ). The component densi¬ 
ties were then obtained using Eq. ( jl4| . In the DM model 
we have a different C = Cj = (co + Ci) ~ 02 , which is in¬ 
dependent of ao- Equation (13 1 was previously obtained 


by Yi et al. Eq] as an improvement over the SMA. The 
breakdown of the single-mode approximation for trapped 
spin -1 condensates in the presence of magnetic field has 
also been theoretically investigated US- 


Provided that ansatz (12 1 holds, distribution (141 can 


be obtained independently from a consideration of C\- 
dependent energy minimization for a ferromagnetic BEC 
as shown in Appendix A. The DM equation is very use¬ 
ful for finding the ground state of a ferromagnetic BEC 
where all density components are non-zero and this pro¬ 
cedure can also be readily generalized to higher-spin cases 
as shown in Sec. |II B| for a spin-2 ferromagnetic BEC. 

For the ground-state of an anti-ferromagnetic BEC 
(ci > 0 ) with non-zero magnetization, energy minimiza¬ 
tion requires that <^o( x ) = 0, viz. Appendix A. Then the 


normalization and magnetization conditions yield 


1 ± TO 

p±idx= —-—, p 0 = 0 . 


(16) 


For to = 0, besides the aforementioned state, there is an¬ 
other degenerate state where all the atoms are in to/ = 0 
component, i.e. p± 1 = 0 and f p gdx = 1 . Unlike in 
a ferromagnetic BEC, ansatz ( fl2| does not hold for an 
anti-ferromagnetic BEC for a non-zero magnetization m 
where different components occupy different spatial ex¬ 
tensions. On the other hand, for to = 0, SMA becomes 
exact in this phase [30], as the ci-dependent term van¬ 
ishes. 

We will derive the analytic model for a ferromagnetic 
BEC using the TFA to the DM equation (131, whereas 


for an anti-ferromagnetic BEC we rely on the TFA to the 
GP equation |2]) with <po(x) = 0 for the same. 


B. Spin-2 BEC 

The dimensionless coupled GP equations for different 
spin components to/ = ±2, ±1 , 0, for a spin-2 BEC can 
be written as [2] 

M±20±2( x ) = di(j)± 2 { x ) + c 0 p(j)± 2 (x) + (c 2 /v / 5)B^2( x ) 

+ ci[F^ ±1 (x)±2F z ^ ±2 (x)], (17) 

M±i</ti( x ) = 7~L(j)±i{x) + cop(/)± i( x ) - (c 2 /V / 5)0^ 1 (x) 
+ Ci [i/3/2F t (/) 0 (x) + F ± (j>± 2 (x) ± F z <j) ±1 {x )], 

(18) 

/i 0 <M x ) = ( x ) +c 0 p4> 0 {x) + (c 2 /v / 5)00 q(x) 

+ c 1 y/3/2[F_<j>_ 1 (x)+F + <t> +1 (x)\, (19) 

where 

F+ =F* = 2(cf)* + 2 ( t ) +i + 4 > -i ( t > -2) 

+ 

F z =2(|</>-t- 2 | 2 — |<?^— 2 I 2 ) + l^+il 2 ~ 10—i| 2 , 

= 20 + 2 0_ 2 - 2(f> +1 (/>_! + 4% 

a/5 

Here the interaction parameters Co = 47rA r (4a 2 + 

3 a4)/(7Zo),Ci = 47rTV(a4 - a 2 )/(7l 0 ),c 2 = 47rY(7o 0 - 
10 a 2 + 3a4)/(7Zo), /r± 2 jM±i; and po are the respective 
chemical potentials. All repeated variables have the same 
meaning as in the spin-1 case. The total density p is again 
normalized to unity. As in the spin-1 case, GP equa¬ 
tions in quasi-2D traps can be obtained by using Eqs. ^ 
and substituting cq = 2A r v / 27r(4a 2 + 3a4)/(7L),ci = 
2TVv / 27r(a4 — a 2 )/(7L),c 2 = 2N\f2n(7ao — 10a 2 + 
3 a4)/(7Z z ) in Eqs. (]I7|-([I9|. Similarly, GP equations in 
quasi-ID traps can be obtained by using Eqs. ( [IT] ) and 
substituting c 0 = 21V(4a 2 + 3a±)lo/(7l 2 z ), C\ = 2N(a± — 
a 2 }lo/(Jll z ),C 2 = 2N(7a 0 - 10a 2 + 3a 4 )Z 0 /(7Z^) in Eqs. 


( 20 ) 

( 21 ) 

( 22 ) 
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In the DM, for a ferromagnetic BEC (ci < 0, C2 > 
20ci) with all non-zero component densities, if we sub¬ 
stitute the ansatz 


= otj(/) dm, j = ±2,±l,0, 


(23) 


in Eqs. (17), (181, and (19), we obtain five independent 
equations for A consistency requirement among 

these five equations for the ci-dependent terms leads to 
the DM equation (131 with C = Cn = (eg + 4ci), pro¬ 
vided that 



(2 ± m) 2 

16 ’ 

(24) 

a±2 = 

u±i \ = 

V4 — to 2 (2 ± m) 

(25) 

8 

Kl = 

y i (4 - m2 )' 

(26) 


with magnetization and normalization conditions 

m= J dx[2(p+ 2 - p- 2 ) + (p+i - p-i)\, (27) 

= 2(|a +2 | 2 -|a_ 2 | 2 ) + (|a +1 | 2 ^|a_ 1 | 2 ), (28) 

1 = |a + 2 | 2 + |a_ 2 | 2 + |a + i|“ + |a_i| 2 + |ao| 2 .(29) 

In the DM model for a spin-2 BEC C = Cn ~ 0,4 is 
independent of the scattering lengths ag and a 2, with 
a 4 playing the role of scattering length in an equivalent 
scalar BEC. With the coefficients a±2, a±i, ag given by 
Eqs. (24l-(p6l the coefficient 0 of Eq. (22) is identically 


equal to 0. The condition 0 = 0 for the ground state 
makes the GP equations © ~( |19[ ) simpler and indepen¬ 
dent of C2- Consequently, the DM equation (13) becomes 


an exact equation for the ground state wave function pro¬ 
vided Eq. (231 holds, e.g. the component wave functions 
are multiples of each other. Our numerical calculations 


show that the condition (231 holds for all magnetization 


to a very high degree of accuracy. 

The coefficients a± 2 , <a±i, ag can also be obtained from 
a minimization of energy along with condition (23), for 


a ferromagnetic ground state with Ci < 0 and C 2 > 20ci. 
An explicit account of the derivation of the coefficients 
a± 2 , a±i, ag, from an energy minimization for a ferro¬ 
magnetic spin-2 BEC is given in Appendix B. 

For an anti-ferromagnetic BEC (02 < 0, C2 < 20ci) for 
any non-zero magnetization m numerical studies show 
that <^±i(x) = 4>g(x) = 0 for the ground state. This can 
also be obtained from energetic consideration as shown in 
Appendix B. The magnetization and normalization con¬ 
ditions (271 and (29) then yield 



FIG. 1: (Color online) Analytic (anal.) and numerical (num.) 
densities of a spin-1 quasi-ID ferromagnetic Sl Rb BEC. The 
number of atoms, scattering lengths and oscillator lengths are, 
respectively, N = 10,000, 00 = 101.8as, a 2 = 100.4as [24j . 
lo = 2.41/rm, l yz = 0.54/xm, here as is Bohr radius. 


Consequently, the magnetization and normalization con¬ 
ditions (271 and (29) lead for these two states 


(i) 


dxp +2 = 


(ii) 


dxp ±2 = 



(31) 

(32) 


For both anti-ferromagnetic and cyclic BECs we will 
derive the analytic models directly from the TFA to the 
GP equations © -(fl9| and not from the DM, whereas 
for a ferromagnetic BEC we will rely on the TFA to the 
DM equation © with C = (cq + 4ci). 


III. ANALYTIC MODEL FOR SPIN-1 BEC 


A. Ferromagnetic BEC 


We derive the analytic model for the ground-state den¬ 
sity of a spin-1 BEC using the TFA to the DM equation 
( fl3| with component densities given by Eq. (141. In the 
TFA the kinetic energy term in Eq. (131 is neglected, 


which is reasonable for a moderate to large positive non¬ 
linear terms, and the BEC density is calculated by equat¬ 
ing the “Hamiltonian” to the chemical potential by 


p = [X 2 /2 + <7pdm], C = Ci = eg + a, 


(33) 


f 2 ±m , , 

Jp ± 2 dx=^—, (30) 

Energy consideration establishes that a cyclic BEC (ci > 
0, C2 > 0) has two degenerate ground states for all 
non-zero magnetization m with (*)^ +1 (x) = <^o( x ) = 
(f >~2 ( x ) = 0, or with (x) = 0, viz. Appendix B. 


thus leading to the TFA density 

Pdm(X) = {ll - X 2 )/ (2C/), X<b =^2p. (34) 

Imposing the condition of normalization f /?dm(A')g?x = 
1, we obtain, in ID, 2D, and 3D, respectively b = 
(3C//2) 1 / 3 , (4C f / tt ) 1 / 4 and (15C//47T) 1 / 5 , provided Cj > 
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FIG. 2: (Color online) Analytic (anal.) and numerical (num.) 
densities of a spin-1 quasi-2D ferromagnetic s 'Rb BEC. The 
number of atoms, scattering lengths and oscillator lengths are, 
respectively, N = 100,000, ao = 101.8 Ob, a ,2 = 100.4os [24] . 
lo = 2.41pm, l z = 0.54pm, here as is Bohr radius. 



FIG. 3: (Color online) Analytic (anal.) and numerical 
(num.) densities of a spin-1 quasi-ID anti-ferromagnetic 23 Na 
BEC. The number of atoms, scattering lengths and oscilla¬ 
tor lengths are, respectively, N = 10,000; oo = 47.36as, 
a 2 = 52.98as [3]; lo = 4.69pm, l yz = 1.05pm. 


0. The component densities are calculated using Eqs. 
(121 and (141. The analytic densities for a quasi-ID spin- 
1 ferromagnetic 87 Rb BEC are shown in Fig. [ljalong with 
the numerical solution of the full coupled GP equations 


§-©• The same for a quasi-2D spin-1 ferromagnetic 
“Tib BEC is shown in Fig. [2] All numerical calcula¬ 


tions are performed using the split-step Crank-Nicolson 
scheme ED with space and time steps 0.025 and 0.00005, 
respectively. 


In the DM model Ci ~ 0,2 plays the same role as the 
scattering length a in a scalar BEC in Eq. 0 . Hence 
in this case the condition of validity of the TFA will be 
Nd2/lo >> 1 . 



B. Anti-ferromagnetic BEC 


In this case, the analytic model is derived by applying 
TFA directly to the GP equation ([2]) with 0o (x) = 0. 
For a non-zero magnetization (0 < m < 1), the irif = +1 
component accommodates more atoms and its spatial ex¬ 
tension 1 -d(+i) is larger than the same of the the rrif = — 1 
component with spatial extension I-d(-i)- Hence for 
It>(+ 1) > x > 1), 0-i (x) = 0 and the coupled GP 

equation (2| for 0±i(x) reduces to a single equation for 
0 + i(x). In the TFA, the kinetic energy terms in the 
GP equation ([2]) are neglected and the densities are cal¬ 
culated by equating the Hamiltonian to the respective 
chemical potentials: 

P±1 = \x 2 /2 + C 0 p\ + Ci(p±i — Pyl), X<l T) (_ 1 ' ) , 

(35) 

M+i = [A r2 /2 + C/p+i] , ^d(+i)>A">Zx>(-i), (36) 


FIG. 4: (Color online) Analytic (anal.) and numerical 
(num.) densities of a spin-1 quasi-2D anti-ferromagnetic 23 Na 
BEC. The number of atoms, scattering lengths and oscilla¬ 
tor lengths are, respectively, N = 100,000; ao = 47.36as, 
a 2 = 52.98as lo = 4.69pm, l z = 1.05pm. 


subject to normalization (161. In the domain 
Eq. (36) has the solution 

P +ipo = l 


;2 _ v2 

L V{+ 1) 


2 C 7 

P+i =lx>(+ 1 )/ 2 - 


> (37) 

(38) 


In the overlap region coupled equations (351 

have the solution 

c 0 (p±i — p T i) + ci(p+i -(- p_i — X 2 ) 

P± il A ) — - 1 -• (o9) 


4c 0 ci 
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The condition p-i(lx>(~i)) = 0 leads to 


P-i = 


Substituting Eqs. 


(co c i)^x>(+i) 2ci lx>(-i) 
2 Ci 


(40) 


and (140)) in Eq. (39I, we obtain 


P- h(A') = 


2co/-p(|i) (^i ^o)^x>(—i) CjX 


4cqCj 


(iL 1 } -x 2 ) 

P-x(X) = 


4 c 0 


X<lx>(- 1). 


X<lx)(- i) 

(41) 

(42) 


The normalization condition (161 for p±i(A) leads to 


k(-i) = ^x>[ c o(l — / C i ] 1 ^ 2+T> \ (43) 

Id(+ l) = Iv[(cq + cim)/Ci] 1 ^ 2+ ' D ' > . (44) 


The densities (371, (411 and B2| with extensions given 


by Eqs. (43)-(44) constitute the analytic model in this 
case. 

These analytic densities for a spin-1 quasi-ID anti¬ 
ferromagnetic 23 Na BEC are shown in Fig. [3] along 
with the numerical solution of the full coupled GP equa¬ 
tions The same for a spin-1 quasi-2D anti¬ 

ferromagnetic 23 Na BEC are shown in Fig. |d] Com¬ 
paring Eq. ([TJ with Eqs. (431-(441, the conditions for 


the validity of TFA in this case are 

N(a 0 + 2a 2 )(l — m) 
Sfo 

-/V[a 0 + 2a 2 + (a 2 - a 0 )m] 
3C 


» T 


» 1, 


(45) 

(46) 


for to/ = —1 and irif = +1 component, respectively. 
The terms on the left side of Eqs. (451 and (46) are the 


measure of the repulsive interactions in to/ = —1 and 
to/ = +1 components, respectively. 

For magnetization to = 0 there is another degener¬ 
ate ground state where all atoms are in the to/ = 0 
component (2). In that case the spin-1 GP equation 


reduces to the DM equation (331 with Cj = Cq and 
p 0 (X) = p DM (X) of Eq. From Eqs. 


the simple criterion for the validity of TFA in this case 
is N(a 0 + 2a 2 )/(3(o) >> 1, which is consistent with the 
fact that for to = 0, Ci term does not contribute to the 
energy of the system. 

The TF analysis shows that the spatial extents of the 
components are different for an anti-ferromagnetic BEC 
for any non-zero magnetization, viz. Eqs. (431-(441, 


which is also manifested by different chemical poten¬ 
tials, viz. Eqs. (381 and (40), whereas these are the 


same for a ferromagnetic BEC. Thus, in the domain 
lv(-i)<X < l-D(+ 1) only component to/ = +1 survives 
for the anti-ferromagnetic BEC effectively separating this 
component from mixed phase in the X < lv(-i) domain. 
The different spatial extents of the components for an 
anti-ferromagnetic BEC also imply that SMA is not valid 


in general except for m = 0, when 1-d(±i) of Eqs. (43)- 
( |44| become equal. The ground states shown in Figs. 3] 
and [4] preserve the symmetry of the trapping potential. 
These symmetric profiles minimize ci-dependent interac¬ 
tion energy m 


E a 


X I [ (P+i ~ P-i) 2 + 2po(p+i + P-i) 
- 4i/p+ip_ipo] dx. 


(47) 


The asymmetric states, where the two phase-separated 
components lie side by side m, do not minimize Ea 
in addition to having more potential-energy contribu¬ 
tion. Hence, they do not emerge as the ground states 
in trapped spinor condensates. The asymmetric states 
can emerge as the ground state in the presence of Zee- 
man energy El El] or spin-orbit coupling ESI ED which 
we do not include in the Hamiltonian. 


IV. ANALYTIC MODEL FOR SPIN-2 BEC 


A. Ferromagnetic BEC 


In this case, the component densities are given by the 
DM equation (231 along with distributions (24)-(26 I and 


the analytic model is derived from the TFATto the DM 
equation (131. Following the procedure of Sec. Ill A for 


a spin-1 ferromagnetic BEC, the TFA densities are again 
given by Eq. ( [34] ) , but now with C = Cjj = Cq + 4ci- 
The component densities are then obtained using Eqs. 
p4|-(|26|). These analytic densities for a quasi-ID spin-2 
ferromagnetic 23 Na and 83 Rb BECs are shown in Figs. 
[5] a) and (b), respectively, together with the numerical 
densities from the full coupled GP equations 

In the DM model C// ~ 04 plays the same role as the 
scattering length a in a scalar BEC in Eq. ([l]). Hence 
in this case the condition of validity of the TFA will be 
Ncl^/Iq >> 1. 


B. Anti-ferromagnetic BEC 


The analytic model here is obtained from the TFA to 
the GP equation ( ]~Tt| ) involving only </>± 2 (x) subject to 
0o(x) = 0±i(x) = 0 with normalization condition (301. 


After neglecting the kinetic energy terms in the GP equa¬ 
tion © the corresponding TFA densities are given by 


P± 2 = X 2 /2 + c 0 p ± 4ci(p +2 - p- 2 ) + 


2 C 2 /?:/ 2 


(48) 


For a non-zero magnetization 0 < to < 2, the to/ = +2 
component extends to a larger domain (X < lx>(+ 2 )) 
than the to/ = —2 component with a smaller extension 
(A < It>(- 2 ),It>(+ 2 ) > ly (- 2 ))- Following the procedure 
presented in Sec. |IIIB| for a spin-1 anti-ferromagnetic 
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FIG. 5: (Color online) (a) and (b) show analytic (anal.) and 
numerical (num.) densities of quasi-ID spin-2 ferromagnetic 
23 Na and 83 Rb BECs, respectively. For 23 Na, the number of 
atoms, scattering lengths and oscillator lengths are, respec¬ 
tively, N = 10,000, ao = 34.9as, a ,2 = 45.8as, a± = 6.45a.s, 
|B] lo = 4.69/rm, l yz = 1.05/im, here as is Bohr radius. The 
experimental value of d 4 (= 64.5cib) has been modified to ac¬ 
cess the ferromagnetic phase of 23 Na (using a Feshbach reso¬ 
nance) from the natural anti-ferromagnetic phase. For 83 Rb, 
the corresponding parameters are N = 10,000, a o = 83.0a_s, 
02 = 82.0 Ob, 04 , = 81.0as, [5] Zq = 2.47/rm, l yz = 0.55/rm. 


BEC, one can calculate lx>(+ 2 ) in ID, 2D and 3D, respec¬ 
tively, as 


, _ [S^l] 1 / 3 , _ [2A] 1 / 4 , [3A] 1 / 5 

1 ( "t 2 ) 20 1 / 3 ’ 2(+2) ~ (57T) 1 / 4 ’ 3(+2) “ (87T) 1 / 5 ’ 

(49) 

where A = (10co + 2c2 + 20cim — cirri '). Similarly, 

i n ID, 2D and 3D, respectively, are 

_ [3 B} 1 / 3 _ [2S] 1 / 4 _ [3S] 1 / 5 

1(_2 ) ~ 20 1 / 3 ’ 2(_2) ~ (57T) 1 / 4 ’ 3(_2) ~ (87T) 1 / 5 ’ 

(50) 


FIG. 6: (Color online) Analytic (anal.) and numerical 
(num.) densities of a spin-2 quasi-ID anti-ferromagnetic 23 Na 
BEC. The number of atoms, scattering lengths and oscil¬ 
lator lengths are, respectively, N = 10,000; ao = 34.9a_e, 
02 = 45.8 ob, 04 = 64.5as [6] ; Zo = 4.69/rm, l yz = 1.05/rm. 


where B = (5co + C2)(2 — m). The analytic TFA densities 
of the spin components irif = ±2 are given by 


P+ 2 PO = 


+ 


20 Cl - A' 2 ) + 2 c 2 (l. 


2 _ / 2 
X>(+2) l V(-2) 


5c 0 ( 


4C//(5c 0 + c 2 ) 


9/2 _ 72 _ v2 

zl X>(+2) l Z>(-2) 


4C//(5c 0 + c 2 ) 


*<b(-2), 

(51) 


P+ 2 (X) = 


X>(+2) 


— X 2 


2 Cn 


’) lv(+2)>X>lxi(-2), 


5 (ll^ - X 2 


P - 2(X) " 4(5c 0 + C2) ’ X -^ ( “ 2) ’ 


(52) 


(53) 


with the extensions l ±2 given by Eqs. (491 and (501. The 


analytic and numerical densities for a quasi-ID spin-2 
anti-ferromagnetic 23 Na BEC are compared in Fig. |6] 
Again, in this case, phase-separated asymmetric profiles 
do not emerge as ground states due to more energy con¬ 
tribution from C 2 -dependent energy term in addition to 
more potential energy as compared to symmetric profiles. 

Comparing Eq. ([!]) with Eqs. (491-( 501, the conditions 
for the validity of TFA in this case are 


-/V[c 0a + 2ci a m + c 2a (2 - m)/10] 
7Zo 

•^[(cpa + c 2a /5)(l - m/2)] 

7L 


» 1 , 


» 1, 


(54) 

(55) 


for rrif = +2 and m/ = —2 components, respectively, 
where c 0a = 4a 2 + 3 a 4 , ci a = 04 — a 2 , and c 2q = 7a 0 - 
10(22 + 3a 4 . 

For m = 0 there is another degenerate ground state 
with the all the atoms in the m/ = 0 component (3]. In 






















this case the GP equation reduces to the DM equation 
(331 with Cn = (co + c 2 /5) and po(X) = pdm(X) of 
(341. A superposition of this solution and the solution 
corresponding to Eq. ( [48] ) with m = 0 will also be a 
degenerate solution. The simpler criterion for the validity 
of TFA in this case is 7V[(7ao + 10a2 + 18a4)/5]/(7Z 0 ) >> 
1, which is consistent with Eqs. (54) and (551 with m = 
0. 


C. Cyclic BEC 


In this case, there are two degenerate ground states for 
all magnetization m with non-zero component densities 
given by Eqs. ( pi) and ( p2) , respectively. The analytic 
models will be obtained in these two cases from the TFA 
to the GP equations 


17j)-( 191. 

The former distribution (311 involves only two non-zero 
components in the GP equation. After neglecting the 
kinetic energy terms in the GP equations the 

TFA densities for the non-zero spin components m/ = +2 
and to/ = — 1 are described by 


P +2 = X 2 /2 + c 0 p + 2ci(2p +2 - p-i), (56) 

p-i= X 2 /2 + c 0 p-Ci{2p +2 ~ p-i). (57) 



FIG. 8: (Color online) Analytic (anal.) and numerical (num.) 
densities of a spin-2 quasi-ID 23 Na cyclic BEC. All parameters 
are the same as in Fig. [7] 


one obtains 

Id(+ 2 ) = Iv[{cq + 2 cim)/Ci] 1 ^ 2+T> \ (58) 

b(-i) = Mc 0 (2 - to)/ 2C i ] 1 ^ 2+v \ (59) 

The normalized densities are given by 


For a non-zero magnetization (2 > to > 0) the to/ = +2 
component has a larger spatial extension (Id<- (- 2 )) than 
the to/ = — 1 component with a smaller spatial extension 
(±7p(_i)> b(+ 2 ) > I'D(-!))■ Following the procedure dis¬ 


cussed for a spin-1 anti-ferromagnetic BEC in Sec. IIIB 



FIG. 7: (Color online) Analytic (anal.) and numerical (num.) 
densities of a spin-2 quasi-ID cyclic 23 Na BEC. The number 
of atoms, scattering lengths and oscillator lengths are, respec¬ 
tively, N = 10,000, ao = 34.9as, a 2 = 22.9as, a 4 = 64.5as 
[Bj, lo = 4.69/tm, l yz = 1.05/Um. The experimental value of 
a 2 (= 45.8as) has been modified to access the cyclic phase 
of 23 Na (using a Feshbach resonance) from its natural anti¬ 
ferromagnetic phase. 


P-1(X) = 


X 2 ) 


3c 0 


X<b(-i), 


(60) 


lv\ ^ C 0^D(+2) 2( c 0 2ci)Z|,(_ 1 ) CuX 2 

P+ 2 (A) = - r - , A <i-D(_!) 


/2 _ \2 
l V{+ 2 ) ^ 

2 On 


6 cqCh 


b(-i)<X<lx>(+2), 


(61) 

(62) 


Equations (60)-(62 I together with extensions given by 
Eqs. (58)-(591 are the analytic densities in this case. 
The analytic and numerical densities for a quasi-ID spin- 
2 cyclic 23 Na BEC are shown in Fig. [ 7 ] For to = 0, SMA 
becomes exact for the cyclic phase of spin-2 condensate, 
as the ci and C 2 -dependent terms in Eqs. (561-( 57) van¬ 
ish. 

Comparing Eq. (jl) and Eqs. (58)-(59), the conditions 
for the validity of TFA in this case are 

IV[4a2 + 3a4 + 2(04 — o 2 )to] 

7 h 

fV[4a2 + 3a4(l — to/2)] 

7L 


»1, 

(63) 

»1, 

(64) 


for to/ = +2 and to/ = — 1 components, respectively. 
For to = 0 the simple criterion for the validity of TFA is 
./V(4a2 + 3ci4)/(7Zo) >> 1, which is consistent with fact 
that only spin-independent non-linearity (co dependent 
term) contributes to the energy of the system. 

Similarly, for the latter distribution (321, after neglect¬ 
ing the kinetic energy terms in the GP equations (|17[)- 
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(19I, the TFA densities are given by 


P ±2 —A 2 /2 + c 0 p ± 4ci(p + 2 — P- 2 ), 
Po =A 2 /2 + c 0 p. 


(65) 

( 66 ) 


This set of equations for p± 2 and p 0 is overcomplete and 
does not determine the densities. However, if we assume, 
consistent with Eq. (|32|, that 


Po 


2(2 + m)p -2 
2 — TO 


(67) 


then we can solve Eqs. (65) for p± 2 and obtain po from 


Eq. (67). For 2 > to > 0 the spatial extent (b(+ 2 )) of 
density 0+ 2 is larger than the spatial extent (lx>(-2)) of 
density p - 2 . Equations (651 can then be be solved to 
obtain 


P+2 (X) - 


4ci(m — 2) (A 2 — S) + cqk(6 + to) 


2p +2 — X 2 


2 On 


64cqCi 


Iv(+ 2 )>X> 1 V (_ 2 ), 


X<lx>(- 2 ) 

( 68 ) 

(69) 


P -2 (X) = 


(to — 2) [c 0 k + 4ci{A' 2 - (5}] 
64cqCi 


X<lx>(- 2 ), 


(70) 


where 5 = p+ 2 + p~ 2 , k = p +2 ~ P- 2 - with the chemical 
potentials p+ 2 and p_ 2 given by 


P+2 


1 ' t >(+ 2 ) (CO 4ci )^X5(_|_2) 4" 8 c 1 ^x>(_ 2) 

~2~ ’ M_2 = 2C// ’ 

(71) 


where lx >(+-2 ) and 1 -d(- 2 ) are the same as lv(+2) and lv(-i ) 
of Eqs. (58)-(591, respectively. After substituting the ex¬ 
pressions for chemical potentials p ± 2 given by Eqs. 
in Eqs. (67 1- (70 1 , we obtain the final densities p± 2 and 
Po as 


P+ 2(X) 


Co 


8Zp( +2 ) - Z|,(_ 2 ) (6 + to) + (-2 + to)A 5 


16coC// 

4ci(2 — m)(Z^, ( _ 2) — A 2 ) 
16coC// 


X<l v{ - 2) , (72) 


P +2 (A) 
Po (A) 


/2 _ y2 

{ X>(+2) ^ 

2 Cu 

2(2 -(- m)p_ 2 

2 — TO 


Z-d(+2)>A>Zx,(_ 2 ), (73) 
X<b(- 2 ), (74) 


P -2 (A) 
Po(A) 


(2 - to) - A 2 ) 

16c 0 

P-2 = 0, Zd(+2)>A ~>l 


, A</p ( _ 2) (75) 

®(-2)- (76) 


The analytic and numerical densities in this case for 
a spin-2 quasi-ID 23 Na cyclic BEC are shown in Fig. [8] 


The criteria for the validity of TFA in this case are again 
given by Eqs. (63)-(64) for to/ = +2 and to/ = —2 
components, receptively! Thus, for a spin-2 23 Na cyclic 
BEC, there are two distinct degenerate ground states as 
are shown in Figs. [7] and [8j The Hamiltonian of the 
spinor BEC is time-reversal invariant, yet the degenerate 
states shown in Figs. [7] and [8] break time reversal sym¬ 
metry. Time-reversal symmetry-breaking states in spinor 
BECs were previously studied m- In cyclic phase too 
the additional potential energy cost rules out the possi¬ 
bility of asymmetric phase-separated profiles as ground 
states. 


V. CONCLUDING REMARKS 


The mean-field GP equation for a spin-1 and spin- 
2 spinor BEC involve three- and five-component com¬ 
plex wave function. Some simplification emerges for the 
ground-state wave function of a spinor BEC. For an anti¬ 
ferromagnetic or cyclic BEC with a non-zero magnetiza¬ 
tion, some of the spin-component wave functions become 
zero, thus reducing the original GP equation with three 
or five components to a system of coupled equations with 
only two or three components, which we call a reduced 
GP equation. For a ferromagnetic BEC with a non-zero 
magnetization the densities of different spin components 
for the ground-state wave function are found to be mul¬ 
tiples of each other. This allows to solve the density 
according to a single GP equation, which we call the 
decoupled-mode (DM) equation, and calculate the densi¬ 
ties of different spin components as multiples of a single 
DM density. These reduced GP and DM equations are 
valid in all spatial dimensions. Here we suggest simple 
analytic models for the ground-state densities of a spinor 
BEC obtained by applying Thomas-Fermi approximation 
to the DM and reduced GP equations. These analytic re¬ 
sults for densities are found to be in good agreement with 
those obtained from the numerical solution of the full 
GP equation for ferromagnetic, anti-ferromagnetic, and 
cyclic spin-1 and spin-2 spinor BECs. Although, we con¬ 
sidered in this paper nearly-overlapping configurations 
of the spinor components, the presence of Zeeman en¬ 
ergy and spin-orbit coupling in the Hamiltonian can lead 
to asymmetric phase-separated configurations nuns] as 
ground states. An investigation leading to the analytic 
densities of the phase-separated solutions would be an 
interesting future work. 
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VI. APPENDIX A 

A. Ferromagnetic spin-1 BEC 

For the ground state of a spin-1 ferromagnetic BEC 
(ci < 0 ), the coefficients ay, can be obtained from a 
minimization of the energy 


(16 1 . If we further assume the DM ansatz (12), the coef¬ 




[ l^l 2 + x 2 p + c 0 p 2 +ci|F | 2 dx, (77) 


j =-1 


where prime denotes x derivative. Assuming that com¬ 
ponent wave functions are given by the DM ansatz (12 1 , 


to minimize energy E we need to maximize the positive 
integral / |F| 2 dx 

J |F| 2 dx = [2\(a* +1 a 0 + a5a_i)| 2 + m 2 ] 1 , (78) 

I = J ^dm (x)dx. (79) 

Now, writing ctj = \otj\e l6 \ we get 

j |F| 2 dx = [2||a + i||a 0 | + |a o ||a-i|e l(e+1+e - 1_20o) | 2 

+m 2 ]l. (80) 

To maximize integral (801, we take exp(0 + i+d_i — 29q) = 
1 and obtain 

j\F\ 2 dx = [2(|a +1 ||a 0 | + |a 0 ||a_i|) 2 + m 2 ]l. (81) 

For a fixed magnetization m and DM function ^>dm> the 
maximization of f |F| 2 dx corresponds to finding the sta¬ 
tionary points of the following “Lagrange” function 

Li(|aj|, Ai, A 2 ) = 2(|a+i||ao| + |ao||a_i|) 2 
+ Ai(l — |cry| 2 ) + A 2 (rn — |a+i| 2 -I- |a_i| 2 ). (82) 


Here Ai and A 2 are Lagrangian multipliers to fix the nor¬ 
malization and magnetization to 1 and m, respectively. 
The stationary points of Li are determined by the fol¬ 
lowing Lagrange equations 

JVo, Sr = 0 ’ Sr = °- ( 83 ) 

o\oij\ a Ai 0 X 2 

with solution together with Ai = 2, A 2 = — 2m. 


B. Anti-ferromagnetic spin-1 BEC 

In case of an anti-ferromagnetic BEC (ci > 0), 
f |F| 2 dx = f(F + F_ + F 2 )dx is minimized by making 
4>o(x) = 0 for any m 7 ^ 0 and the densities satisfy Eq. 


ficients aj can be obtained from a minimization of (801. 
For this, we take exp(0 +1 + 9_i — 2 9 0 ) = —1 and obtain 

J\F\ 2 dx = [2(|a +1 ||a 0 Mao||a-i|) 2 + m 2 ]l. (84) 

Following the procedure discussed for a ferromagnetic 
BEC, one can minimize / |F| 2 dx under the twin con¬ 
straints of fixed norm and magnetization and, in agree¬ 
ment with Eq. (16), obtain 


1 1 ± m 

|Q(-j-i| \l — , Q!q 0* 


(85) 


VII. APPENDIX B 

A. Ferromagnetic spin-2 BEC 

For a spin-2 ferromagnetic BEC (ci < 0 ,C 2 > 20ci), 
the energy is given by 


r 2 

/ [ + x2 p + C 0P 2 + Cl IF1 2 + C2 1 ©I 

J A — O 


f _ N 

I L ^ 

i=-2 


dx. 


For a ferromagnetic BEC, the energy minimization cor¬ 
responds to a maximization of the Ci-dependent term 
f |F| 2 dx. We find that this automatically minimizes 
the C 2 -dependent term f |@| 2 dx to zero. Assuming the 
DM ansatz (23) we seek the coefficients ay which maxi¬ 
mize f |F| 2 dx. Following the procedure for ferromagnetic 
spin-1 BEC, we can write 

/|F| 2 dx = 2{|a + 2 ||a + i| + |a_ 2 ||a_i| 

— 6-1 — 81 +S 2 ) | _|_ y/gg*(0O — 201+02) 

{l a +i|l a o| + |au||a_i|e l( - e - 1 ~ 2 eo+e, b} + m 2 l. ( 86 ) 


To maximize (861 we take all exponential factors in this 


equation to be +1 and obtain 


IFrdx = 


+ a _ 2 a_J 


+ \/6(|a+i||a 0 | + |a 0 ||a_i|)| +m 2 l. (87) 

For a fixed m (2 > m > 0) and I, the maximization of 
/ |F| 2 dx corresponds to finding the stationary points of 
the following Lagrange function 

L 2 (|ay|,Ai, A 2 ) = 2(|a +2 ||a + i| + |a_ 2 ||a_ 1 |) + \/6|a 0 | 

(l^+iI + l a -il) + Ai(l - \ a j\ 2 ) + A 2 

3 

(m - 2|a +2 | 2 - |a + i| 2 + |a_i| 2 + 2|a_ 2 | 2 ). 
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Here Ai and A 2 have the same meaning as in Eq. (821. 
The stationary point which maximizes f |F| 1 2 <ix thus 
yields Eqs. (241-(26 I together with 


Ai = - 




m 2 ) 


Ao — — 


m 


^(4- TO 2 ) 


( 88 ) 


Using Eqs. (24)-(261, we find that these ay's also min¬ 
imize f |0| 2 uxto 0 which guarantees that the state so 
obtained is the ground state. 


B. Anti-ferromagnetic spin-2 BEC 


The stationary point, which maximizes j |0| 2 dx, yields 


\ot±2\ = 

Ai = 


\/2 ± 


TO 


a±i = ao = 0, (91) 

— TO. 


^/(4 — TO 2 ) ’ 


A 2 — 


2 ^ (4 — to 2 ) 


■ (92) 


Using Eqs. (911-(921, we find that f |F| 2 dx has the min¬ 
imum value rnrX which guarantees that the state so ob¬ 
tained is the ground state. 


Similarly in anti-ferromagnetic subdomain (02 < 0 and 
C 2 < 20ci), the energy minimization corresponds to a 
maximization of the C 2 -dependent term J|0| 2 dx. As¬ 
suming DM ansatz (161 and aj = \oij | exp (iOj) we have 


J |©| 2 dx = \2\a +2 \\a- 2 \ ~ 2\a +1 \\a- 1 \e l ^ +e - 1 - 62 - e - 2 

\2gi(20o —@2 — 0-2) \ 2 j^ 




(89) 


To maximize integral (891 we take the first exponential to 


be —1 and the second exponential to be +1. For a fixed 
to (2 > to > 0) and I, the maximization of f |0| 2 dx cor¬ 
responds to finding the stationary points of the following 
Lagrange function 

Lg(\aj\, Xi, X 2 ) = (2|cr_|_21|cr—21 + 2|a+i||a_i| + |«o|") 


C. Cyclic spin-2 BEC 

For a spin-2 cyclic BEC ci > 0 and C 2 > 0, energy 
minimization involves minimization of both f |F| 2 dx and 
f |0| 2 dx to their respective minimum values m 2 X and 0, 
respectively. From equations (861 and (|89|), one can see 


that for 0 < to < 2, consistent with Eqs. (|31[|-(32 1, there 
are only two possibilities for the ground states: 


Wl«±2| = 


2 ± TO 


K*o = 


\/4 — TO 2 


4 7 8 9 10 11 y/8 ’ 

provided that exp(20 o — 9+2 — 9 - 2 ) = —1, and 


a±i = 0 , 


+Ai(l - ^2 \ a j\ 2 ) + - 2|a +2 | 2 (i)|a+ 2 | = 


\/l + m , , V2 - m _ 

, \Qt— 1| , Q^±l Q^o 0- 


-|a +1 | 2 + |a_ 1 | 2 +2|a_ 2 | 2 ). (90) 
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